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We investigate the properties of subharmonic gap structure (SGS) in superconducting quantum 
contacts with normal-electron resonances. We find two distinct new features of the SGS in resonant 
junctions which distinguish them from non-resonant point contacts: (i) The odd-order structures on 
the current-voltage characteristics of resonant junctions are strongly enhanced and have pronounced 
peaks, while the even-order structures are suppressed, in the case of a normal electron resonance 
being close to the Fermi level, (ii) Tremendous current peaks develop at eV = ±2Eo where Eq 
indicates a distance of the resonance to the Fermi level. These properties are determined by the 
effect of narrowing of the resonance during multiple Andreev reflections and by overlap of electron 
and hole resonances. 
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I. INTRODUCTION 

Resonant tunneling often plays an important role in 
current transport in transmissive superconducting junc- 
tions. The presence of impurities in the ituanel barri- 
ers may strongly affect Josephson tunnelinJi □ as well as 
Andreev transport under applied voltageoH due to the 
resonant tunneling through localized impurity levels. A 
mechanism of resonant tunneling was assumed to be re- 
sponsible for the unusual properties of junctions with 
disordered semiconducting barriersQ and grain bound- 
ary Josephson junctions of high Tc materials!! Further- 
more, in clean superconductor -semiconductor junctions, 
mobile electrons confined between the Schottky barriers 
form resonant states which determine specific properties 
of such junctionso'E2l. Similar resonant states are also im- 
portant in recently fabricated ballistic superconductor- 
2D electron gas (2DEG) structures^ where they are 
formed by the electron reflections by the gate potential. 
Moreover, the possibility to create quantum point con- 
tacts and quantum dots in such structures allows one 
to investigate quantumjesonant transport through well 
resolved resonant levelsE3. Quantum resonant transport 
was also observed in metallic datsL3, and in contacts con- 
taining single wall nanotubesll-jll^ . 

The properties of dc Josephson current in quantum 
resonant junctions as well as Andreev ow^Htization have 
been discussed in various publications^ - !!*]. The sub- 
harmonic gap structure (SGS) in resonant junctions and 
the effect of the resonance on multiple Andreev reflec- 
tions (MAR) is less investigated. Meanwhile- detailed 
theory of SGS in quantum point contactsE2Tc3 in com- 
bination , wj | th precise experiments on controllable break 



junctions!: 



was found to be a powerful tool for in- 



vestigation of intrinsic properties, of the atomic-size con- 
tacts. Our preliminary resultsEZI have shown that SGS 
in resonant junction drastically differs from the SGS in 
non-resonant junction; similar results were obtained by 
different methods in Ref.E3. In this paper, we present 
a detailed study of interplay of MAR with Breit-Wigner 
resonances in quantum junctions. 

The structure of the paper is the following. In Sec- 
tion II we derive equations for the inelastic scattering 
amplitudes in resonant junctions. In Section III we 
discuss properties of the normal electron resonance in 
the proximity region between the superconducting elec- 
trodes. Results of numerical calculations of the current- 
voltage characteristic for different positions and widths 
of the resonance are presented in Section IV. Section V 
is devoted to perturbative analysis of resonant SGS. 



II. SCATTERING AMPLITUDES. 

We will consider a junction consisting of a ballistic 
normally conducting region separated from the super- 
conducting electrodes by tunnel barriers, Fig. 1. The 



length of the junction L is assumed to be smaller than 
the coherence length, and therefore the distance between 
normal resonances will exceed the superconducting en- 
ergy gap, vf/L > A [vf is the normal electron Fermi 
velocity, h — 1). We will also assume that the reso- 
nances are well separated, T <C vp/L, where T is the 
resonance half-width, and that Coulomb charging effects 
do not dominate in the subgap voltage region, Ec < 2A, 
where Ec is a Coulomb gap. 

We Y^lVl a PPly Landauer-Biittiker scatter 
ing theorynjEy extended to superconducting junctionsEll 
for calculating the current. In voltage biased supercon- 
ductive junctions, the quasiparticle scattering is inelastic 
due to time dependence of the superconducting phase 
difference across the junction, <fi(t) = 2eVt, and the scat- 
tering state wave functions consist of linear combinations 
of harmonics (sidebands) with energies E n = E + neV 
shifted by integer number of quanta eV with respect to 
the energy E of the incoming wave. Below we will con- 
sider one transport mode in the junction and choose the 
scattering state wave functions in the left (L) and right 
(R) electrodes having the form 
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where are (non-normalized) two-component elemen- 
tary solutions of the Bogoliubov-de Gennes equations, 



u£ = 



V2\ v n e 



3 ±7„/2 



(2) 



In this equation 

^ \E n \+L 



ia n VA 2 - El \E n \ < A 



a n = sgn(-E„), fc„ = ^2m(E F ± a n £„). 

Index j = 1 — 4 in Eqs. (|l|) labels scattering states of 
the electron- and hole-like quasiparticles incoming from 
the left (j — 1, 2) and right (j = 3, 4). The form of wave 
functions in Eqs. (|l]), (Q) assumes that superconducting 
electrodes serve as equilibrium quasiparticle reservoirs, 
and that the potential difference between the reservoirs 
is absorbed into the time-dependent factor e 1 ^^ in Eq. 
( |lb| ) due to appropriate choice of the gauge. 

To match the wave functions in Eqs. (m) we will ap- 
ply a transfer matrix technique. In the present case of 
an inelastic scattering problem, the connection between 
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ipL and is non-local in time, and the corresponding 
transfer matrix TfL-. mixes the sidebands, 
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The matrix T„ m is a 4 x 4 matrix defined on a space 



of wave function coefficients, A 
(B + ,B-), 
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The transfer matrix in Eq. (|3|) can.be expressed, similarly 
to the case of unbiased junctionstj, through a transfer 
matrix T(E) associated with elastic electron scattering 
by the normal junction. Let us introduce auxiliary nor- 
mal regions between the superconductors and the junc- 
tion, the length of which normal regions will be put equal 
to zero at the end of the calculations. The wave function 
in the normal region has the form 
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where — y/2m(Ep ± E n ) is the normal electron 

wave vector. The wave functions Eq. (|^) in the left and 
right normal regions are related as 
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We note that the transfer matrix T(E) describes scat- 
tering of the normal electrons by the actual potential of 
the junctions at a given voltage, i.e. it includes effects of 
potential deformation due to applied voltage, T(E; V). 

Continuous matching at the left SN interface yields in 
quasiclassical approximation k n « « fcf, 
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A matching condition at the right NS interface is derived 
in a similar way but an additional time-dependent factor 
gicr^evt j n ^b) must be taken into account. The 

latter gives different equations for upper (electron) and 
lower (hole) components of the coefficient vectors: 
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In this equation, are projectors on upper/lower vector 
components. 



Combination of Eqs. (^J)-(||) gives the following equa- 
tion for the transfer matrix in Eq. ([j|) 
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Normal electron transfer matrix enters this equation with 
different arguments ±E n . This energy difference intro- 
duces effects of electron-hole dephasing during quasipar- 
ticle propagation through the junction. In non-resonant 
short constrictions, the energy dispersion of the trans- 
fer matrix is negligible, and equation (OL is equivalent 
to the matching equation derived in RefEJ. In resonant 
junctions (and also in long SNS and SIS junctional) de- 
phasing effects are important. 

The matching equations ([}]), and (^|) can be written in 
an equivalent form, 
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Applied to the scattering state wave functions in Eqs. 
(|l), it yields the following recurrences for the scattering 
amplitudes: 

e^ 2 M&)+e--- /'(«)„ = 
T(E n ) [e-^ 2 5 (n+1)0 (|:) +e-^ +1 / 2 (/) n+i 
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Analytical solutions of the recurrences in Eqs. ([Il]) can 
be presented in chain-fraction form (see A pmci jidix A) 
similar to the case of non-resonant junctionsEEil. 



III. MODEL FOR RESONANCE. 

Now we will specify the transfer matrix for the reso- 
nant junction. We will restrict ourselves to symmetric 
junctions, T\\ = T 2 * 2 — 1/d and = T* 2 = r/d and as- 
sume the Breit-Wigner resonance form for transmission 
and reflection amplitudes d and r respectively, 
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The position of the resonance level E r as well as the 
resonance half-width T are generally dependent on the 
applied voltage. However, while the subharmonic gap 
structure is affected in an essential way by the position 
of the resonance, the dependence on the resonance width 
is less important. Thus we will assume T — const. We 
will not specify the voltage dependence of resonance level 
position, but rather present the current as a function of 
two variables: driving voltage and resonance position, 
I(V, E r ). The current voltage characteristics can then be 
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reconstructed from such a dependence by specifying the 
E r (V) dependence determined by self-consistent distri- 
bution of the electric potential across the junction. 

The normal electron resonance, being confined between 
superconducting electrodes, possesses specific properties 
which will be important for further analysis of the res- 
onant MAR. Since the transfer matrix T(E) enters the 
recurrences for the scattering amplitudes in Eq. ( [ll| ) at 
two different energies ±E, the proximity resonance con- 
sists of two, electron and hole, resonances situated sym- 
metrically with respect to the Fermi level, E = ±E r . 
Within the adopted approach, the current is calculated 
by using the scattering amplitudes defined in the super- 
conducting electrodes [see further Eq. (|l5|)] and the recur- 
rences in Eq. (|Ti"| ) are formulated for these amplitudes. 
Although equivalent, such an approach is different from 
the discussion of MAR amplitudes in the normal region 
of the junction (see, e.g.E3). Within our approach, the 
non-superconducting region of the junction is considered 
as a black box and is represented by the transfer ma- 
trix T(E). Due to the different choices of gauge in the 
left and right electrodes, the resonance is seen from the 
left and right electrodes at different energies [cf. Eqs. 
H and (|)]. Indeed, the resonances are seen from the 
left electrode at E = ±E r , i.e. quasiparticles incoming 
from the left undergo resonant transition if E = ±E r , 
while the resonances are seen from the right electrode at 
E = ±(E r + eV), as shown in Fig. 2. In the scattering 
diagram in Fig. 2c, the resonance therefore is presented 
with two segments: E n = E r <-> E n+ \ = E r + eV for the 
electron resonance, and E n — —E r «-> E n -\ = — E r — eV 
for the hole resonance. 

There is a symmetry between the scattering states orig- 
inating from the left and right electrodes: 
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with an analogous relation for the second pair of scat- 
tering amplitudes. In Eq. (jl^), E = E r + eV/2 is the 
distance of the normal resonance level with respect to 
the midpoint between the chemical potentials in the left 
and right electrodes. Equation ([l3]) leads to a symme- 
try property of the current which is an even function of 
the resonance position Eq: I(V,Eq) = I(V, —Eq). Below 
we will indicate the resonance position by means of the 
energy Eq and abbreviate the Breit-Wigner amplitudes 
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IV. DC CURRENT. 



In the quasiclg 
the current reads 
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The current in Eq. dig ) is calculated using transmitted 
states (in the right and left electrodes for scattering states 
j = 1,2 and j = 3,4 respectively), and it consists of 
contributions from all odd sidebands. By virtue of the 
symmetry equations 
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directly following from Eqs. (O) (analogous relations 
hold for the scattering states j — 3, 4) and symmetry 
equations (|l^), the current in Eq. ( |l5| ) can be expressed 
through the sideband contributions 
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of one single scattering state {j = 1, index j is omitted), 
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where K n = K n (-£ n , -V). 

Equation Eq. ( p^8[ ) together with the recurrences in 
Eqs. (|ll| ) provide a basis for numerical calculation of the 
current. The calculation of scattering amplitudes should 
obey the boundary condition at ± infinity where the am- 
plitudes approach zero. The simplest way to obtain such 
solutions is to iterate the recurrences from large \E n \ to- 
wards E. The correct solution will then grow exponen- 
tially and numerically "kill" the solution growing at in- 
finity. By this procedure one gets the correct scattering 
states for each incoming quasiparticle at every energy. 

The results of numerical calculation of current-voltage 
characteristics (IVC) are presented in Fig. 3 for differ- 
ent values of resonance level position Eq = const. This 
particular case corresponds to a perfectly symmetric dis- 
tribution of the electric potential across the junction with 
Eq indicating the departure of the resonance level from 
the Fermi level in equilibrium (V = 0). The IVC with 
the resonance level situated at the Fermi level, i^o = 0, 
shows an onset of the single-particle current at eV = 2A 
accompanied by a current peak caused by large density of 
states near the superconducting gap [see below Eqs. ( pij - 
Pq)]. Such behaviour of the single-particle currept has 
been observed in the experiments on metallic dotald and 
carbon-nanotube junctions^. A striking feature of this 
IVC is the absence of current structure at eV — A, while 
the structure at eV = 2A/3 is pronounced, consisting of 
a peak similar to the structure of the single-particle cur- 
rent. Calculation of the IVC at lower voltage, presented 



4 



in Fig. 4, shows the same feature - only odd subharmonic 
gap structures are present. 

If the resonant level departs from the Fermi level and 
Eq 7^ 0, the single-particle current onset shifts towards 
larger voltage, eV > 2A, and the current peak broadens. 
A striking feature in this case is the development of a 
huge current peak at voltages lower than the position of 
the structure of single-particle current. This peak, as- 
sociated with resonant pair current (see below Sec. V), 
appears as soon as E > A/2 and is situated at voltage 
eV = 2E which coincides with position of the resonant 
current onset in normal junctions. If the resonance level 
departs far from the Fermi level, Eq A, the IVCs in 
the subgap voltage region eV < 2A approach the form 
typical for non-resonant point contacts, as could be ex- 
pected, while a broadening of the resonance, r » A, 
gives rise to SNS-type IVCs, as shown in Fig. 5. 

A complete description of the current in resonant junc- 
tions is given by the function I(V,E ), as already men- 
tioned in Section III. A plot of this function is presented 
in Fig. 6. The IVCs plotted in Figs. 3-5 correspond to 
horizontal cuts (E = const) of the plot in Fig. 6. In Fig. 
6a, the light wedge-like region at eV > 2A corresponds 
to the resonance single-particle current. The resonant 
peak of the pair current is seen as the light streaks di- 
rected along the lines Eq = ±eV/2, the structure starts 
at eV = A. Fig. 6b presents a similar plot for the re- 
gion of small voltage, eV < A. The picture shows quite 
complex structure of the current consisting of wedge-like 
plateaux of the resonant current as well as of light streaks 
corresponding to current peaks. 

In order to interpret the features of the IVCs one needs 
to analyze the properties of the side-band currents K n 
presented in Eqs. (|l7|) and (A7). 



V. DISCUSSION 

A convenient expression for analysis of the subhar- 
monic gap structure is derived in Appendix B: 



are hidden in Eq. Jig)]. This is done by proper renor- 
malization of the sideband currents K n — > K n presented 
in Appendix B, the equation for K n being given in Eq. 
flB9p. We will now develop a perturbative analysis of 
the current in the limit of small width of the resonance, 
r <C A, and zero temperature. 



A. Single-particle current 

The single-particle current is given by the first term in 
Eq. (|l9|). In accordance with Eq. (BE), it has explicit 
form 
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Pi is defined in Eq. ( JBlOj ). This current has no contri- 
bution from Andreev reflections and it has only one res- 
onance. It is sufficient to consider only scattering states 
incomin g fr om the left [the first term in curly brackets 
in Eq. (p0|)], the resonance equation in this case being 
E^ = 0. The resonance is only involved if it belongs to 
the integration interval. This determines the resonance 
region eV/2 > A + \E \ in the plane (V, E ) (region / 
in Fig. 7). The resonant scattering diagram is shown in 
Fig. 8a. 

In non-resonant junctions, the currents K^ n in Eq. 
@ have the singularities which are responsible for the 
current onset at eV — 2A and subharmonic gap struc- 
ture. In resonant junctions, the singularities are washed 
out due to strong electron-hole dephasing and the res- 
onant transmissivity is simultaneously renormalized. In 
the case of a single-particle current in Eq. (|2(]), the onset 
of non-resonant current is caused by zeros of the function 
Pi. Calculation o f Pi for the resonant junctions by using 
the rule in Eq. (B1C) and retaining only the resonant 
scattering amplitude yields 
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In Eq. ( |T9| ) only contributions of processes of creation 
of real excitations (transitions across the gap E > A — > 
E n < — A) are retained mhich are responsible for the sub- 
harmonic gap structures while a contribution of ther- 
mal excitations is omitted. Furthermore, the sum over 
the sideband currents in Eq. (18) is now rearranged in 



order to explicitly separate the contributions of all in- 
elastic channels [contributions of even inelastic channels 



A 4p2 
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where r„ = r|_E„|/£ n . Equation (|2l| ) shows the transfor- 
mation of the resonant tunneling probability in the super- 
conducting junctions: the resonance width is broadened 
due to superconducting density of states E/£. Taking 
into account Eq. ( pl| ) and similar equations for the other 
terms in Eq. (|20|), we present the single-particle current 
on the form of the Landaucr formula, 
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with the effective single-particle transmission coefficient, 
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Similar equation have been derived in Ref.Ea using a dif- 
ferent method. 

Equations ( |22| ) and (|2^) determine the current in the 
wedge region in Fig. 6a. In the limit of T — > 0, the 
resonant current reads 
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V- ^{eV+f - 4A 2 + V+^eV-f 
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where eV± — eV ± 2 1 -Eq | - This equation quantitatively 
describes the single-particle current feature in Fig. 3. 
The current has maximum at the wedge edges and de- 
creases at large eV approaching the value for the resonant 
current in the normal junction In = eT (see Fig. 3), 

, / 2|£ ° I+A ,eV = 2 { A + \E \) 

h=lN \ vmm+A) (25) 
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The current peak is the result of enhancement of the 
effective width of the resonance in Eq. (|23| ) at low en- 
ergy C = Equation ( p5| ) is everywhere applicable ex- 
cept of the wedge vertex, E — 0, eV — 2A, where the 
current grows without limit. In fact, the current con- 
sists of the peak and turns to zero at eV = 2A due to 
shrinking of interval of integration in Eq. ( p2| ) in this 
region. The maximum current is achieved when the inte- 
gration interval becomes comparable with the resonance 
width, eV - 2 A ~ T v / A/{eV-2A). These arguments 
yield estimate for the maximum current at eV — 2A, 

(h)max ~ /AKA/r) 1 ^. 



B. Pair current 



The pair current has form 
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Restricting again the consideration with quasiparticles 
incoming from the left, we find that this current is con- 
tributed by two resonances, Eq = and E\ = 0, which 
simultaneously enter the integration interval within the 
region Il\ in Fig. 7 (region II2 corresponds to resonant 
quasiparticles incoming from the right). Therefore, the 
resonant pair current only appears if the normal reso- 
nance is sufficiently far from the Fermi level, Eq > A/2, 
while at i?o < A/2 the current is non-resonant within the 
voltage interval A < eV < 2A. This means in particular 
that the onset of the pair current at eV = A is small: 
I2 ~ /Ar(r/A) 3 if Eq — 0. In regions 77, the pair cur- 
rent undergoes resonant enhancement, I2 ~ I/v(r/A) 2 
due to independent contributions of two separate reso- 
nances (Fig. 8b), each contribution being described by 
the equations similar to Eqs. (E2h, (23). 



The most interesting phenomenon in the resonant pair 
current is the overlap of the resonances occurring along 
the lines eV = ±2|£o| in Fig. 7. The overlap of the 
resonances produces a huge current peak near these lines 
(seen as light streaks in the Fig 6a; we note that these 
lines correspond to the position of the onset of resonant 
current in normal junctions). The scattering dia gram for 
this case is presented in Fig. 8c. Applying Eq. ( |B1C| ) for 
calculation of P2 and retaining both the resonant ampli- 



tudes d and d_ 2 , we obtain 
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Substituting Eq. (|27]) into Eq. (|26|) for the current and 
collecting the contributions of all scattering modes, we 
find 
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Equation (|29|) shows a remarkable similarity to the res- 
onant transmissivity of Schrodinger three-barrier struc- 
tures: the probability to leak outside the superconduct- 
ing gap through the sidebands n = and n — — 2 (Fig. 
8c) corresponding to probability of tunneling through 
side barriers, while the probability of Andreev reflection 
by the sideband n = — 1 corresponding to transmissivity 
of a central barrier. Such three-barrier structures have 
been investigated, e.g. in connection with normal . elec- 
tron transport properties of coupled quantum dotsoEl 
The strong overlap of the resonances is explained by the 
fact that shift of the resonances is proportional to T 2 
due to Andreev reflection, according to Eq. (p9|), while 
the resonance width is proportional to the first power 
of r (the quantity T_i is equal to zero at the lines 
eV = ±2\E \). 

In the vicinity of the lines eV — ±2|75o| an d in the 
limit r — > 0, the pair current has following form 
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r 2 eVy/(evy 



(eV-2\E \) 2 [(eVy 



A 2 ]+r 2 [2(eV) 2 -A 2 ]' 
(30) 



(We notice that this formula is valid at all voltages eV > 
A because the side band n = — 1 is inside the energy gap 
if eV ~ ±2|i?o|). Equation ( |30| ) describes the current 
peak in Fig. 3, the height of the peak 
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being comparable to the magnitude of the resonant 
single-particle current, in particular, (^Jmai = In /2 at 
E > A. 

According to Eq. (|3^) the resonant pair current tends 
to zero at large voltage eV 3> A,Eo, which means that, 
rigorously speaking, there is no resonant excess current. 
However, if the resonance is far beyond the gap, \Eq\ 3> 
A, the current may strongly deviate from the current in 
the normal junction in the region A <C eV <C 2|£'o| be- 
cause the single-particle current is non-resonant in this 
region, while the pair current is resonant. Such an effect 
is particularly pronounced in the junctions where the res- 
onance level follows the chemical potential of one of the 
electrodes, E (eV) ± eV/2 « e = const. The IVC in 
this case corresponds to cuts in the plot in Fig 6a paral- 
lel to the light streaks. In such a case, the peak of the 
pair current is very broad, and even transforms into a 
plateau with a sharp onset at eV = A (e = 0), as shown 
in the inset in Fig. 10. The magnitude of the cur rent at 
the plateau can be found directly from Eq. (A7) when 
assuming Eq = e ± eV/2 and eV — oo, 
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D t = r2 /[((£ T e) 2 + T 2 ]. This current as function of e 
is shown in Fig. 10. 

There is an interesting difference between the property 
of the resonance in the single-particle current and that of 
individual resonances of the pair current. To be specific, 
let us consider the resonance Eq: in the pair current 
this resonance is more narrow because the quantity r_i 
is imaginary and causes a resonant shift rather than a 
contribution to the resonance width. The physical reason 
for this squeezing of the resonance is that direct leakage of 
a quasiparticle through the side band n = — 1 is blocked, 
and the only escape from the resonant region into the 
continuum is through the states of the side band n = 0. 



C. High-order currents 



respectively. The side resonances at Eq , EZ 2 = are 
characterized by an effective transmissivity similar to 
the effective transmissivity of the resonances of pair cur- 
rent (times additional factor ~ T 2 ). The contribution of 
these resonances in the current is therefore estimated as 
Is ~ 7jv(r/A) 4 . The central resonance E^ 2 = is much 
more narrow. Indeed, in this case (Fig. 8d), direct leak- 
age of the resonant particle into continuum is blocked at 
the both sidebands n = —1,-2, and the particle can es- 
cape only through the sideband states n = 0, —3, travers- 
ing the junction one more time. The central resonance 
determines the current in the vicinity of the threshold 
eV = 2A/3, E = 0. 

Calcu latio n of the quantity P3 in region III 2 according 
to Eq. (|B1C| ) yields 



DqD± 2 D~_ 2 



A 4 fnf_ 



\&&- Z EE- 



E^ 



r 



(34) 



where E+ 2 = E+ 2 + i(r_i + r_ 2 )/2 + 0(T 2 ) and T = 
r ^ -A 2 /4^| 2 , f_ 3 = r_ 3 £l 2 A 2 /4|£- 2 | 2 . According 
to Eq. (pj), the resonance width is of the order of 
r ~ r 3 which yields giant enhancement of the current 
I3 ~ iAr(r/A) 2 , exceeding by two orders of T the con- 
tribution of the side resonances. Such narrowing of the 
central resonance occurs in the quadrangle region in Fig. 
7 bounded by the edges of the resonance region III 2 and 
regions II. The current in this region has a form similar 
to the one in equation (E2h, 



3eV — A 



dE D 3 (E), 



A 



with the effective resonant transmissivity 



D 3 (E) 



3r r 



E+ 2 -i r + r_ 3 /2 



(35) 



(36) 



In the limit of small r —> 0, the current reads 



/3 = 6 e r f_ 3 / r + f_ 



E=\E \+3eV/2 



(37) 



The effect of the resonance narrowing is even more 
important for the third order current, 



3eV-A 
A 



-I' A /"" - a dE \E-M-3 
A 3 



e 7 v?_3 



DqD1 2 DZ 2 



Pa 



P 5 



The third order current has three resonances at 
Eq , E^ 2 , EZ 2 = which belong to the interval of in- 
tegration within the regions III\, III 2 , III3 in Fig. 7, 



The phenomenon of resonance narrowing provides the 
explanation for the absence of current structure at volt- 
age eV = A, namely the dominance of the third order 
current 23 at the threshold of the pair current. The cur- 
rent in Eq. ([57]) is responsible for the light wedge-like 
region at eV < A in Fig. 6b. Similarly to the case 
of single-particle current, the third-order current in Eq. 
d37| ) has a peak at the edges of the wedge with the height 
+ (£^0 — ► — Eo) f • (33) increasing proportionally to (eV — 2A/3) -1 / 2 towards 
the vertex of the wedge, eV = 2A/3, E = 0. This 
growth is again limited due to interplay between shrink- 
ing integration interval and growing resonance width, 
eV - 2 A/3 - T 3 [A/(eV ~ 2A/3)] 1 / 2 . This estimate 
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gives a height {I^) ma x ~ In{^/ A) of the current peak 
at eV = 2 A/3. As one may see in Fig. 6b, there are 
no current structures at the edges eV = 2 (A — |-Eb|) of 
the above mentioned quadrangle where the the narrow 
resonance of three-particle current dies: this is because 
of the resonant pair current emerges at the same lines, 
giving rise to a gradual cross over between three-particle 
current and pair current both having the magnitude of 
the order of I N (T/A) 2 . 

The phenomenon of resonance narrowing results in 
enhancement of central resonances in all higher odd- 
order currents, giving rise to the current peaks at 
eV = 2A/(2fc + 1), E = with the height I max ~ 
J/v(r/A) 2fc ~ 1 . The magnitude of the current between 
neighbouring peaks is I ~ iAr(T/A) 2fc . Also, the overlap 
of narrow resonances of the even-order currents near the 
lines eV — ±2\Eq\ yields current peaks with the height 
I ~ /Ar(r/A) 2fe within the intervals 2A/(2fc + 2) < eV < 
2 A /2k. These current peaks are clearly seen in Fig. 6b 
in the form of light streaks. 



VI. CONCLUSION. 

In conclusion, we have considered effect of the nor- 
mal electron resonant tunneling on the subharmonic gap 
structure (SGS) in mesoscopic superconducting junc- 
tions. In non-resonant tunnel junctions, the SGS con- 
sists of sharp onsets and narrow peaks of the current 
at voltage eV — 2A/n. In resonant junctions, SGS is 
considerably modified depending on the position of res- 
onance level with respect to the chemical potentials of 
the electrodes. If the resonance level is situated exactly 
in the middle between the chemical potentials of elec- 
trodes, the odd-n current structures are tremendously 
enhanced while the even-n current structures are not af- 
fected by the resonance. This enhancement is explained 
by narrowing of the resonance during multiple Andreev 
reflections. When the resonance departs from the mid- 
point between the chemical potentials of electrodes, new 
current structures appear at eV = ±2Eq in a form of cur- 
rent peaks. This feature results from overlap of electron 
and hole resonances. 

In our calculations, the Coulomb charging energy was 
assumed to be smaller than the superconducting gap, and 
the charging reflects were neglected. fai©K»criments on 
metallic dotsEJ and carbon nanotubesEJO'Ed, the oppo- 
site situation has been observed with the Coulomb charg- 
ing energy exceeding the superconducting gap which led 
to suppression of the subgap current. The charging en- 
ergy in quantum transport experiments can be reduced 
by enhancing capacitance of the resonant structure, e.g. 
by using substrates with large dielectric constants. This 
will allow direct application of our results to such struc- 
tures. Another way would be to use high-T c materials 
for fabrication of superconducting electrodes for the nan- 
otube experiments. Our theory is applicable to ballistic 



plane junctions with large capacitance such as resonant 
junctions in high mobility S-2DEG-S devices and atomic 
plane junctions in layered cuprates (intrinsic Josephson 
junctionsEZl) . Current-voltage characteristics of such mul- 
timode junctions can be obtained on the basis of our 
theory by summation of contributions of all transport 
modes. 
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APPENDIX A: DERIVATION OF CURRENT 



Following the method of RefBJ'Eil, we eliminate the 
Andreev scattering amplitudes a n and b n from Eqs. ( |Tl] ) 
for the scattering state j = 1 and obtain a closed set of 
equations for the normal amplitudes c„ and /„, 

c 2n + V 2n+1 f2n+l + V 2 ^_ 1 f 2n ~l = -r^ S nO (Al) 

Ann 



/2n+l + V^ +2 C2n+2 + V 2n C2n — 0. 

The coefficients in these equations are 



v 2n — a 2n e /"2n+l, 



v. 



2)1 



v. 



2n+l 



-a 2 „a 2 „_ 1 d 2 -*e-^"+^"- 1 )/ 2 /S 2 
4„*e (72 " +72 "+ l)/2 /S 2 „, 



a -(72»+72»-l)/2 /n 



l^2n, 



where the quantities H n are defined as 



Ho , — r + * P T2»-i 
—2n-l — '2ri-2 e 



-*„-72r, 



* e — 72n-l 



(A2) 



In non-resonant junctions, the functions H„ approach 
r£„ since energy dispersion of the reflection amplitude 
is negligibly small. SGS in non-respnant junctions is 
caused by zeros of the functions and renormaliza- 

tion of these functions in the resonance case, — » H„, 
is the reason of considerable difference of the SGS in 
resonant junction compared to the SGS in non-resonant 
junctions. We solve equations (Al) by introducing ratios 
S2n = Ca n /f2n-i and S 2n +i = hn+x/c2n and expressing 
/„ through c , 



fn — Y\_ &i c 0) 
i=0 



(A3) 
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Co be ing related to S±i by virtue of the first equation in 
Eq. (Al). By introducing chain-fractions Z ni 



o7+7i 



Z a = l- «*) 2 ^— - (dg*) 2 „ - , (A4) 



-I2n— 72n-l 



S2n^2ri-1^2n-l 

g72n — 1+T2n — 2 



"2n-l"2ri-2^2n-2 



Zln—1 ' 



one can rewrite Eq. (A3) for f n , e.g. with negative 
sideband index n = — 2fc — 1<0, in the following form 



2/c+l 



(A5) 



An equation for 6-2fe-i follows from Eqs. (|ll| ) and ( A5), 

(Z3 + g7-2fe-2 \ 
1 + + ' 2 i'" 2 7 /-afc-i, 
r -2/c-2"-2fc-2^-2fc-2 / 

(A6) 



where £± = |d±| 2 . 

Collecting the norm al and An dreev transmission am- 
plitudes in Eqs. (|A5| ) and (A6) and substituting them 
into Eq. (03) for the current, we finally get 



K n = -|/„| 2 cosh(i?e7„) 



1 - e 



2Re~/ r , 



(i-Ati) 




The corresponding equation for positive n > is obtained 
from Eq. ( A7) via the substitutions 7 — * —7, D + — > 
and n — 1 — * n + 1 . 



APPENDIX B: TRANSFORMATION OF if* 



Equation ( A7 ) is not cpavenient for analysis of the sub- 
harmonic gap structureEil. 



The current structures are 
caused by the processes of creation of real excitations 
during across-the-gap transitions E — > E n < —A. The 
currents of other side-bands (E^-> A) perfectly cancel 
each other at zero temperatureEJ. Furthermore, there 
is rigorous balance between Andreev and normal channel 
currenta_among the states lying within the superconduct- 
ing gapcj, which provides successive drops of the current 
with decreasing applied voltage. All of these features 
are not explicitly seen in the currents K n in Eq. (A7). 
Moreover, the current structures related to the creation 
of real excitations in even-order sidebands arc hidden in 



Eq. ( |l7| ) which consists of contributions from only odd 
sidebands. 

To overcome this diffi culty , we will transform the side- 
band currents-dn Eq. (A7) following the method sug- 
gested in Ref.Eil. Assuming the transparency Dn-i m Eq- 
(A7) to be small, we find that the current K n to lead- 
ing order is proportional to sinh(i?e 7Tl ) ~ 0(E% — A 2 ). 
This observation will allow us later [in Eq. (19)] to sepa- 
rate the contributions from states below, E n < —A, and 
above, E n > —A the energy gap. Having separated out 
the leading term, we rewrite Eq. (A7) in the form 



K n = \fn\ 2 cosh(i?e 7 „ 
( 2sinh(i?e 7 „)e fle7 ' 1 - 



a2Re~f n 



|S n _iZ„_i 



(Bl) 



F n = \~ n Z n \ 2 - 2Re (e<r+3 n Z n ) - D+e 2Re ^. (B2) 



Equation (B2) possesses a similar property: it is propor- 
tional to 0(Ef l — A 2 ) to leading order with respect to D. 
Sepa rati ng out this leading term, we further transform 
Eq. (p3) into the equation 



F„ = -2sinh(2i?e 7lI )-^^- „G n - U (B3) 

pn-i Ai-i| 

where 

Gn-i = \^n-iZ n -i\ 2 + 

+ 2Re(e-<-ir-E n - 1 Z n ^ - D~e- 2Re ^-\ (B4) 

One more transformation, 



G n -x = 2sinh(2i?e7„_i) 



3n-2-Zn-2 



-F n . 



(B5) 



accomplishes the cycle, yielding the quantity F in Eq. 
(B2) with shifted index. Performing repeate dly such 
transformations, we get for the current in Eq. (A7) the 
following expansion 



= 0(El - 


-A 2 )Q n 


+ 




M(K-i 


-A 2 )e" 




cosh(i?e 7n )Q„_i + 




-A 2 )e" 




+2Re 7n - 1 cosh ( i?e7n )Q n _ 2 + 




-A 2 )e" 


-Rej r 


+2fe 7 „_ 1 -2Ke 7 „_ 2 cosh ( i?e7n )Q n 



(B6) 



In this equation, the quantity Q n is defined as 



A 4 P„ 



'k-1 



Do [I[D± 2i DZ 2i )D± 2k 



1, \n\ = 2k 
DZ 2k , |n|=2fc+l 



(B7) 



where 
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H 

Pn = l[ \Z-rZ- 



(B8) 



i=0 



Collecting together all terms with similar ^-functions, 
we can finally rearrange the sum in Eq. (|lS|): 

oc 

12 Kn £ = ^ - A 2 )Q„^„, (B9) 

n—odd n—1 

where is given by the recurrence equation 

<p n -i = 1 + exp[(-l)" +1 2i?e(7„)]^„, <p_ x = 1. 

Far from resonance, the quantity P n may cause strong 
singularity in the sideband current due to the presence 
of zeros in the functions 3 n , and accounting for the fac- 
tors Z n is .absolutely necessary for regularization of the 
singularity^. In the resonant case, the functions S n do 
not turn to zero because of strong electron-hole dephas- 
ing, f$_£ r~ , and the quantities Z n can be omitted from 
Eq. (p38|) in the limit of narrow resonance V <§C A, 



Pn 



(B10) 



The role of Z n in this limit reduces to cancellation of 
the terms in the product ( B10| ) which are proportional 
to the squared resonance a mplit udes (d^) 2 ; this is de- 
noted by the prime in Eq. ( B10| ). The presence of the 
resonant denominators in equation (A.2) for S gives rise 
to renormalization of the normal electron transmission 
coefficients in Eq. (|l9|) for the current. 
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A. 




-A 



FIG. 1. Scheme of the double-barrier junction. Upper part: 
dark regions - tunnel barriers; shadowed region - normal con- 
ductor. Lower part: energy diagram showing one normal res- 
onant level inside the superconducting gap. 



\Ef*reV 



■E r +eV 




eV/A 



FIG. 3. IVC of symmetric resonant junctions, a) IVC of 
normal junction (dashed line) and superconducting junction 
(solid line), Eq/ A= 0, 1.2 (left and right curves respectively), 
T/A = 0.05. b) IVC of superconducting junction, Eq/A = 0. 
(1), 0.6 (2), 0.9 (3), 1.2 (4). 



a) 



b) 



r E+eV 

/, — \ 

) E-eV 

c) 

FIG. 2. Energy diagram of resonant junctions under ap- 
plied voltage, a) Resonance in normal junction; the distance 
to the chemical potentials of the left and right electrodes is E r 
and E r + eV respectively, b) Electron and hole resonances in 
superconducting junction; the resonance offset is different in 
the left and right electrodes with respect to the global chem- 
ical potential after equalizing the chemical potentials of the 
electrodes by means of the gauge transformation, c) Resonant 
transition in the scattering diagram. 



Q-4 
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eV/A 



FIG. 4. Subharmonic gap structure on IVC of symmetric 
resonant junction, Eo — 0, T/A = 0.2. The structures appear 
only at eV = 2A/n with odd n. 
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eV/A 

FIG. 5. IVC of symmetric resonant junctions with Eq — 
and r/A SE {0.1,0.4,0.7,1.0,2.0,3.0,10.0} (from bottom to 
top). The current is normalized by Io = eA/hn. 




eV/A 



Log(l/l N ) 



FIG. 6. Intensity plots showing the dependence of the cur- 
rent I(V, Eq) on the applied bias voltage V and the resonance 
position Eq. (a) Single-particle current and pair current at 
eV > A (r = 0.2A). (b) Pair current and high- order cur- 
rents at eV < 2A (r = 0.05 A). 
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2A 



-eV/2 
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FIG. 7. . Resonant regions in the plane (eV,Eo,) for side- 
band currents I n . I, II and III are resonant regions for sin- 
gle-particle current, pair current and three-particle current 
respectively. The resonant regions are bound by dotted lines 
E = ±(A - neV/2) (labeled with ±n). Bold dashed lines 
show positions of double resonances. 
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-A 
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a) 



b) 




0- 



/ 



c) d) 

FIG. 8. Resonant scattering diagrams; bold lines show 
resonant transitions, a) Resonance in single-particle current; 
b) single resonance in the pair current; c) double resonance 
in the pair current; d) central resonance in the three-particle 
current. 
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FIG. 9. Spectral density of the pair current as a function 
of bias voltage for T/A = 0.05. The current peak appears at 
the crossing point of two single resonances. 



1.5 




e/A 



FIG. 10. Current for asymmetric junction at large volt- 
age (eV — 1000A as a function of the position of the 
resonance e = \Eo\ — eV/2 for different resonance width 
r/A G {0.01, 0.21, 0.41, 0.61, 0.81, 1.01} (from top to bottom 
at the left side of the figure). Inset shows IVC for e = 0: the 
resonance level coincides with the chemical potential of one 
of the electrodes. 
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